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Abstract 

It is shown that curved and flat helical double twisted liquid crystal (DTLC) in blue phase, 
can be unstable (stable) depending of the sign, negative (positive) of sectional curvature, 
depending on the pitch of the helix of the nematic crystal. In both cases Cartan tor- 
sion is presented. It is also shown that the instability or stability depends on the value 
of the pitch of the helix in nematic crystals. Frank energy stability A similar result us- 
ing the method of Frank energy stability in the twist of cholesteric liquid crystal was given 
by Kiselev and Sluckin [PRE 71(2005)], where the free twist number determines the equi- 
librium value of the cholesteric liquid crystals (CLC) pitch of the helix. As a final ex- 
ample we solve the geodesic equations in twisted nematics with variable pitch helix and 
non-constant torsion. This non-Riemannian geometrical approach, also seems to unify two 
recent analysis of cylindrical columns given by Santangelo et al [PRL 99,(2007)] and the 
curved crystal endowed with torsion, given by Vitelli et al [Proc Nat Acad Sci (2006)]. 
Stability of toroidal curved surfaces were also previously considered by Bowick et al [PRE 
69,(2004)], as an example of curvature-induced defect. Investigation of the Lagrangean in- 
stability may be useful in the investigation of HIV viruses and proteins. PACS numbers: 

02.40.Hw:Riemannian geometries 
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I Introduction 



Following the steps of success of the applications of the differential geometry of curves and 
surfaces [1], through Riennannian geometry in Einstein theory of gravitation [2], as well as 
the success of the use of Cartan non-Riemannian geometry endowed with a torsion tensor, 
to explain alternatives theory of gravitation such as Einstein-Cartan gravity [3], many mathe- 
matical and crystal physicists, have applied these geometries to the explanation of disclination 
(Riemann curvature) or dislocation (non-Riemannian Cartan torsion), defects in crystals [4]. 
Even recently Riemannian geometry has been applied in the investigation of twisted magnetic 
flux tubes in plasmas helical flows MHD with and without vortices [5, 6, 7]. More recently 
these paths, two recent papers on the stability of geodesies in the curved and flat cylindrical 
liquid crystals in curved surface with defects [8] and columnar phases [9] have appeared in the 
literature. In the first paper Vitelli et al [8] have discussed the instability driven by a purely 
geometric potential in Gaussian curvature of the stress law covariant generalization. In the 
second Santangelo et al [9] have been considered the instability or stability of geodesies in 
the geometry of columnar phases of nematic liquid crystal, in the analogy with the focusing 
process of a lensing ray in optics or gravitational lensing. Recently we also use the sectional 
developed by Kambe [10] curvature method to investigate the Lagrangean instability of a Cou- 
ette type flow, viscous and sheared. Actually Kambe [10] has previously investigate Couette 
planar flows stability, and showed that Riemann curvature tensor would exist in the case of 
existence of pressure, even a constant one. The stability of the Couette planar flow is also 
obtained by him, using the symmetry in the case of symmetric scalar stress potential and 
by making use of the technique of Riemann sectional curvature, where the negative sectional 
curvature indicates instability of the flow, in the Lagrangean sense. In this brief report this idea 
can be easily transported to investigate the instability of DTLC nematic blue phase, where 
as shown by Pansu and Dubois-Violette [11] the presence of torsion is taken from granted 
from the very beginning. As a last example we modify the constant torsion twisted nematic 
geometry with variable pitch and compute the geodesies for this CLC. Therefore we investigate 
the Lagrange instability in two simple examples, in analogy with the works of Vitelli et al and 
Santangelo et al, namely of a Riemann-flat torsioned DTLC , while in the second example 
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a curved crystal is considered. The usefulness of the investigation of the Riemannian and 
non-Riennannian geonnetries in physical allows one also building models for analog gravity as 
the case of non-Riemannian geometry of vortex acoustic flows [12]. Other curvature- induced 
defect in toroidal Riemann geometry were also investigated by Bowick et al [13]. Yet another 
example of the use of Gaussian curvature on continuum elastic model topological defect in 
Riemannian manifolds was investigated by Giomi and Bowick [14]. The paper is organized as 
follows: Section II presents a brief review of sectional Riemannian curvature method in the 
coordinate-free language which helps to unify previously examples of instability. Section III 
presents the examples of curved and flat DTLC with totally skew torsion Lagrange stability and 
the dependence of the instability on the value of the pitch of helix, a similar result obtained by 
Kiselev and Sluckin [15]. Section IV presents the computation of geodesies of twisted nematic 
geometry. In section V conclusions are presented. 

II Sectional Riemannian curvature Instability 

In this section, before we add we make a brief review of the differential geometry of surfaces 
in coordinate- free language. The Riemann curvature is defined by 



where XeTM is the vector representation which is defined on the tangent space TM. to the 
manifold M.. Here represents the covariant derivative given by 



which for the physicists is intuitive, since we are saying that we are performing derivative along 
the X direction. The expression [X, Y] represents the commutator, which on a vector basis 
frame ei in this tangent sub-manifold defined by 



R{X, Y)Z := VxVyZ - VyVxZ - V^x,y]Z 



(II.l) 



VxY = (xv)y 



(II.2) 



X — Xkik 



(II.3) 



or in the dual basis dk 



X = X^di 



'k 



(II.4) 
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can be expressed as 

[X,Y]^{X,Y)'dk (II.5) 
In this same coordinate basis now we are able to write the curvature expression (II. 1) as 

R{X, Y)Z := [R'jkpZ'X''YP]di (II.6) 

where the Einstein summation convention of tensor calculus is used. The expression R{X, Y)Y 
which we shall compute bellow is called Ricci curvature. The sectional curvature which is very 
useful in future computations is defined by 

K^--{X, Y) (11.7) 

where S{X,Y) is defined by 

s{x,Y) := < x,y >2 (ii.s) 

where the symbol <, > implies internal product. In the non-Riemannian (NR) case, the torsion 
two-form T(X, Y) is defined by 

T(X, Y) := h;SIxY - VyX - [X, Y]] (II.9) 

where V is the non-Riemannian connection [7] endowed with torsion. As in EC theory [5] 
the geodesic equation does not depend on torsion; only Jacobi deviation equation depends 
on torsion which is enough for investigate the role of torsion on stability. Since the Jacobi 
equation is given by 

^ = [||V,-e.|r-ir^^(t,e,)||J|| (11.10) 

where \ \ej\ \ — 1 and J is the Jacobi field, representing the separation between geodesics,while 
t is the geodesic tangent vector. Here K^^(X, Y) is given by 

K^^{X, Y) = K^'^'^{X, r) + 2 < T(X, Y)SyX > (11.11) 

Here as we shall see bellow the geodesic equation is VyF = simplified this expression. Note 
from this expression that the instability, or separation of the geodesies in the flow 

g > (11.12) 
implies that K^'^^"^ < which is the condition for Lagrange instabillity. 



Ill stability of Curved and fiat Liquid Crystal with 



torsion 

In this section let us address the examples of Riemann-flat TDLC and the curved LC. In the 
first case we consider the CLC, where twist of molecules exists in just one direction, which a 
ordered phase is given by the director field 

= (III.13) 

— cospx (III. 14) 

— sinpx (III. 15) 

where q :— ^ \s the wave number and p is the pitch of the helix. In the blue phase local 
constraint of minimum energy is double twist , which physically means that twist occurs in 
more than one direction in space. In the Euclidean frame this equation is 

5kn^ + pskipu'S^P ^ (III.16) 

In the Riemann-flat manifold, where the Riemann tensor R{X, Y)Z vanishes, or in components 

R'jki ■■= (III. 17) 

In the frame where the covariant derivative of the basis (cj) vanishes or Vg^ej — 0, where 
from expression (11.9) the Cartan torsion reads 

T{ei,ej)^-[ei,ej] (III.18) 

The double twist Levi-Civita connection is 

Vifcn^' = dkfi^ + Pifeff = (III. 19) 

Now computing the general Riemann-Cartan curvature tensor 

R^'^iCi, Cj) = Ve-,Ve-, - Ve-,Ve- - V[S„e,] (111.20) 

Using the expression for the Cartan torsion 2-form T{X, Y) in terms of the Lie bracket above, 
the RC curvature reduces to 

R^'^iei, ej)ej = -V[e- ,e-,]e,- = ([e,, e,-].V)e,- (III.21) 
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therefore to compute this expression one computes the Lie bracket in terms of Cartan torsion 
as 

T{ei, Cj) = (ei.V)ej - (e}.V)ei = [cj, Cj] = diCj - djCi = r%-4 (III.22) 
substitution of these components of torsion into expression (111.21) yields 

i?^^(e-;, e.y = ([el, e,].V)e^ = (r\4.V)e'^ = (T^-^.e') = T^.T'^^el (III.23) 

Thus the cholesteric blue phase of liquid crystals possesses the sectional curvature 

K{ei, ej) =< R{ei, ej)e^, >= T'^'^Tkij = -pe^jke''" = -6p' < (III.24) 

since this curvature is negative, the instability of the CLC is proved in this case. In a more 
general Riemannian manifold we shall prove now that depending on the value of the pitch 
the CLC can be stable in the Lagrangean sense. This can be done as follows. First we shall 
compute the Riemann tensor with the condition of geodesic VyF = is given, which simplifies 
the curvature expression to 

R{X, Y)Y = -Vy[VxY + pYAX] + p{XAYV)Y (III.25) 

where now, we have used the covariant derivative 

VxF := VxY + pYAX] (111.26) 

To simplifies computations we assume that the term Vx^ vanishes which after some algebra 
leads to expression 

R{X, Y)Y = [-p^ + p]YA{XAY) (III.27) 
Finally from this expression we can write the expression for the sectional curvature as 

K{X,Y) =< R{X,Y)Y,X >= [-p'^+p] < YA(XAY),X >= [-p'^ + p]a\\X\\^ (III.28) 

Thus this time the sectional derivative can be positive or the twisted nematic crystal can be 
geodesically stable, depending on the value of the helix pitch. 
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IV Geodesies of twisted nematic with variable pitch 
hehx 



Let us built the Riemann metric in analogy to Katanaev and Volovich [16] as the deformation 
of the crystals = ?T.(j_i) where this deformation represents the perturbation on the metric 
Qij given by 

9ij = Sij + (IV.29) 

where S^j is the delta Kronecker symbol and the vector where {i = 1,2,3) is the director 
field of nematic liquid crystal, we consider here that Riemann curvature tensor vanishes as 
in Einstein theory of absolute parallelism, which implies that the anti-symmetric connection 
(torsion tensor) is given in terms of the metric by 

Tijk — gjk,i — gji,k (IV.30) 

Here the comma denotes partial derivatives with respect to the lower index. Substitution of 
expression (IV.29) into (IV.30) yields the following torsion vector 

Tk = [V^nfe - Skidivn)] (IV.31) 

therefore and the Weitzenbock condition for teleparallelism on the curvature Riemann tensor 
Rijki(r) — .where F is the Riemann-Cartan affine connection. Later on we shall make an 
application of this formula to a special case of nematic crystal. In the meantime let us com- 
pute the geodesies equations for the corresponding metric of the liquid crystal. The geodesies 
equations 

+ Tii^x^x^ = (IV.32) 

and 

X' + 5'^nk,ijX^x^ = (IV.33) 

where r*^. = |<^*'[^ij,fe + 9ik,j — gjk,i\ and we use the Euclidean 3D metric Sij to raise and 
lower indices. Let us now apply these ideas to the pure twist geometry of the nematic liquid 
crystals.where the director field is now given by the components 

n, = cose(y) (IV.34) 
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and 

= sin9{y) (IV.35) 

Where 9 is the twist angle and the planar crystal is orthogonal to the y-coordinate direction. 
Substituting expressions (IV. 33) and (IV. 34) into (IV.35) one obtains the following components 
of the torsion vector T^- — 

Ty{e) = -dyidivn) (IV.36) 

and 

TM = -V^n, (IV.37) 

Since divn — nx,x + n'y,y + n:^,z = , there is no torsion vector component along the orthogonal 
direction. Thus the only non-vanishing torsion component of torsion reads 

T,{e) = -[cosd{ — f + sinOy^] (IV.38) 

Since local equilibrium conditions on the Liquid crystals yields [11] 

dO 

— = constant = K (IV.39) 
dy 

Thus equation (IV.38) reduces to 

T,{e) = -K^cose (IV.40) 

This last expression tell us that the twisted geometry of the crystal leads to a 3D helical 
torsion. Geodesies of the twisted geometry of the liquid crystal leads to the results 

x + nlyy' = Q (IV.41) 

and 

y = (IV.42) 

and finally 

z + n^yf^O (IV.43) 

Here the dots means derivation with respect to time coordinate. A simple algebraic manipu- 
lation yields 

X + Ki^sind{y) = (IV.44) 



and 

z + K^^cose{y) = (IV.45) 

To simplify matters let us solve these equations in the approximation of small twist angles 
9 «< where above equations reduces to 

X + Ki^e = (IV.46) 

and 

z + Ki^^O (IV.47) 
Therefore from expression (IV.47) one obtains the following solution 

z = -KiH + c (WAS) 

where c is an integration constant. To solve the remaining equations we need an explicit form 
of the twist angle with respect to time ,this can be obtained from the integration of the 
expression ^ — which yields 

e^Ky + f (IV.49) 
where / is another integration constant. Substitution of y = Kot + d into this equation yields 

e^mt + g (IV.50) 

where f,m and g are new integration constants. Thus one obtains 

X = -{K\^mt + Kig) (IV.51) 

Integration of this expression yields 

x{t) = ~{at^ + + 'yt + 6) (IV.52) 
substitution of the expression t = into (IV.52) yields 

X = ~iay' + P'y^ + j'y + S) (IV.53) 

where all the greek letters represent new integration constants. Therefore the trajectory of the 
test particles is represented by a third order polinomial curve. Trajectories of test particles in 
domain walls are in general parabolic curves. After we finish this letter we hear that Dubois- 
Violette and Pansu [11] have considered a similar application of teleparallelism to cholesteric 
Blue Phase of liquid crystals. Nevertheless in their paper Cartan torsion is constant 
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V Conclusions 



One of the most physical question in liquid crystal physics is the investigation of the stability of 
flows in the Euclidean manifold using the sectional curvature of geodesic deviation. In this 
paper we discuss and present the contribution of Cartan torsion tensor and its role in curved 
and Riemann flat liquid crystals in the blue phase. The Lagrange instabilities in several CLC 
cases are computed and it is shown that instability may depend upon the value of twisting 
number or the pitch of the helix of nematics. The models discussed here may also be useful 
in building the analog models of stability of geodesies in non-Riemannian theories of gravity. 
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